Atomic Interactions in Precision Interferometry Using Bose-Einstein Condensates 



O 

(N 

< 
(N 



C/2 

-t— > 



I 

c 

o 
o 



> 

cn 

o 



X 



Alan 0. Jamison, 1, J. Nathan Kutz, 1,2 and Subhadeep Gupta 1 

1 University of Washington Department of Physics, Seattle, Washington 98195, USA 
2 University of Washington Department of Applied Math, Seattle, Washington 98195, USA 

(Dated: January 14, 2013) 

We present theoretical tools for predicting and reducing the effects of atomic interactions in 
Bose-Einstein condensate (BEC) interferometry experiments. To address mean-field shifts during 
free propagation, we derive a robust scaling solution that reduces the three-dimensional Gross- 
Pitaevskii equation to a set of three simple differential equations valid for any interaction strength. 
To model the other common components of a BEC interferometer — condensate splitting, manip- 
ulation, and recombination — we generalize the slowly-varying envelope reduction, providing both 
analytic handles and dramatically improved simulations. Applying these tools to a BEC interfer- 
ometer to measure the fine structure constant, a (Gupta, et al., 2002), we find agreement with 
the results of the original experiment and demonstrate that atomic interactions do not preclude 
measurement to better than part-per-billion accuracy, even for atomic species with relatively large 
scattering lengths. These tools help make BEC interferometry a viable choice for a broad class of 
precision measurements. 



I. INTRODUCTION 

Interferometry with atomic beams or laser-cooled gases 
enables a number of high-precision measurements for ap- 
plied and fundamental research [l[: rotation sensing 
measuring atomic properties , and testing foundational 
principles of general relativity [4J and quantum mechan- 
ics Q , to name a few. Bose-Einstein condensates (BECs) 
have properties valuable to precision experiments: long 
coherence lengths (and correspondingly narrow momen- 
tum distributions) and high density. These "atom lasers" 
thus have tremendous potential for high-accuracy mea- 
surements. 

Several promising experimental results have been 
shown with free-space [6j-[8| and confined (in traps or 
wave guides) [9Hll| BEC interferometers. While much 
theoretical work has been done to understand the effects 
of atomic interactions in confined interferometers [l2T - 
[l5j . the quality of the trapping potential tends to domi- 
nate considerations for high-accuracy measurements. In 
this work we present theoretical techniques to deal with 
atomic interactions in free-space BEC interferometers 
where these interactions are the fundamental hurdle for 
high- accuracy measurements. 

In this paper we derive and apply two complementary 
techniques which together address the general problem 
of interaction effects: Robust scaling solutions correctly 
model BEC free expansion at any interaction strength, 
and a decomposition method efficiently models multi- 
path interferometers. Together these techniques allow 
modeling of a broad class of condensate interferometry 
experiments, both analytically and with simplified simu- 
lations. With these techniques it is now possible to ex- 
plore a wide range of parameters in a short time to opti- 
mize experimental design as well as to run high-precision 



lamisonaiauw 



edu 



simulations of chosen parameters on commodity PCs. 
Our results provide key new tools for precision BEC in- 
terferometry. 

We consider a large class of experiments in which a 
condensate is split (i.e., placed in a superposition of mul- 
tiple momentum states), subjected to some differential 
phase shift, and then recombined. To keep the discus- 
sion focused, we will apply these techniques to a BEC 
interferometer to measure the fine structure constant, a 
[fjj, which uses standing- wave laser diffraction pulses as 
mirrors and beam splitters for the condensate [16]. Our 
techniques agree both with full numerical simulations and 
the results of the first-generation contrast interferometry 
experiment. 

For the class of experiments considered in this pa- 
per, collective excitations such as solitons, vortices, and 
Shockwaves [l7j are confined to the short time windows 
of condensate splitting and recombination and can be 
monitored with numerical simulations fl8j ] . We find no 
such excitations in any of the simulations reported here. 
Further, densities and atom numbers are low enough to 
keep elastic collisions infrequent and below the thresh- 
old for bosonic stimulation 19]. This leaves the mean- 
field shift — the energy shift of an atom due to its inter- 
actions with all other atoms in the BEC as the most 
important interaction effect for our work. While mean- 
field effects may, in principle, be reduced using Feshbach 
resonances 20] , doing so introduces systematic Zeeman 
shifts and the need to make nonzero magnetic fields ultra- 
stable and uniform. Further, the residual interaction ef- 
fects need analysis, for which our techniques are valuable. 

The paper is organized as follows. In section[TT]we dis- 
cuss a robust analytic approximation for calculating the 
mean- field shift. In section UTT1 we develop a generalized 
slowly- varying envelope approximation. In section llVl we 
apply our tools to identify experimentally accessible pa- 
rameters such that mean-field effects do not preclude a 
sub-part-per-bllion (sub-ppb) measurement of u; rcc and 
a. 
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II. ROBUST SCALING SOLUTION TO THE 
GROSS-PITAEVSKII EQUATION 



In discussing atomic interaction effects, we will con- 
fine ourselves to the zero-temperature Gross-Pitaevskii 
description of a BEC, which is quite accurate for 
many interferometry experiments [171 ] . This approx- 
imation assumes that all atoms are in the conden- 
sate and that the many-body wavefunction has a form 
i>{xi , . . . , Xpft) = Ylj <fi(%j j t) ■ The Gross-Pitaevskii 
equation (GPE) is the mean-field approximation for <p 
and its temporal dynamics: 
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where g — Airh 2 a s (-/V at — l)/m, a s is the s-wave scattering 
length and iV at is the total atom number. 

In this section we describe an approximate solution to 
the GPE which is accurate for both small and large in- 
teraction strength. Our scaling solution builds from a 
previous result valid in the Thomas-Fermi (TF) approx- 
imation, derived by Castin and Dum [2l|. The TF ap- 
proximation describes condensates where the interaction 
energy is much larger than the kinetic energy and may 
be considered as a semi-classical limit. Castin and Dum 
derived a "dynamical TF" approximation by explicitly 
considering this semi-classical nature, using an analogy 
to the expansion of a classical interacting gas to find a 
self-similar scaling solution to the GPE. Using an alter- 
nate derivation we find scaling solutions accurate near the 
origin for any interaction strength, which reduce to the 
result of Castin and Dum in the TF regime. In particular, 
our scaling solution reproduces the exact solutions both 
in the g — > oo (TF) limit and for g = (noninteracting). 

For precision interferometry experiments, bringing the 
mean-field shift as low as possible without compromis- 
ing the advantages of a BEC is essential, so the residual 
interaction effects will likely be well outside the regime 
where TF is applicable. As will be seen below, interac- 
tion effects for such BECs can still be large enough to 
spoil the results of precision interferometers. 

We will consider the time evolution of a condensate in 
free space [22j]. Suppose the initial state has some given 
form (j)Q (x). Now, consider an ansatz of the form 



</> (x, t) = 



t)> A 2 (t) ' A 3 (t) 
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i8(x,t) 



(2) 



where Xi, 9, and 4>q are all real valued. In the follow- 
ing, we will use the notation yt = xt/Xi to simplify the 
appearance of equations. This gives d/dxi — > X~ l d/dyi 
and means there is some time dependence to <po(y) due to 
the Ai's. We substitute this ansatz into (p} and separate 
into two equations — one for the imaginary part and one 
for the real part. The equation for the imaginary part 
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Setting each of the two expressions in parentheses equal 
to zero we find the following condition on 9: 



3=1 



(4) 



With this form for 9, the real part of the GPE yields 
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Thus, if ([5]) can be satisfied exactly for a given <f>o, the 
ansatz gives us an exact solution to the GPE. 

We first look at the two extremes of interaction 
strength. For the TF limit of a condensate released from 
a harmonic trap at t = 0, we have 
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(g-too), (6) 



where f2j/(27r) is the trapping frequency in the j direc- 
tion and \x is the chemical potential in the trap. The 
second term on the righthand side of ([5]) becomes negli- 
gible, recovering the original result of Castin and Dum. 
In this limit, <j5j> can be solved exactly. As noted above, 
this result may be interpreted as the expansion of a clas- 
sical, interacting gas driven by pressure due to repulsive 
interactions. From this perspective, the second term can 
be viewed as the first quantum correction, correspond- 
ing to an extra quantum "pressure" due to wavepacket 
dispersion. 

In the opposite limit of a noninteracting gas (or a single 
atom) released from a harmonic trap, (f>o is gaussian and 
again ([5]) can be solved exactly. Thus, we can expect ([2]) 
subject to (|4]) and §5§ to give good results for very small 
or very large g. To extend this solution to all g values, 
we view the above equations as an expansion in y. Thus 
we find solutions valid for all interaction strengths in a 
region near the center of the condensate. 

The signal in most interferometry experiments will be 
dominated by the highest density region when the arms 
overlap. This is true because the majority of atoms will 
be found here. Also, with tiny offsets from perfect over- 
lap, as are inevitable in experiments, the rapid phase 
oscillations in the wings will largely cancel any contri- 
butions far from the slowly- varying peak density region. 
Furthermore, since the primary goal is to understand how 
to suppress the phase shift due to mean-field effects, set- 
ting an upper bound is key. Accurately knowing the 
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FIG. 1. (Color online) Accuracy of scaling solutions at center of BEC. The phase at the center of the condensate is plotted 
versus expansion time for (a) high density, (b) intermediate density, and (c)low density. The red circles show results from full 
numerical simulations, while the solid blue curves show our scaling solutions using the initial condensate wavefunction calculated 
in the trap. For reference, the TF scaling solutions are plotted as dashed black curves and the noninteracting solution is plotted 
as thick green curves. Neither the TF nor the noninteracting solution is sufficient for high-accuracy measurements at any of 
these densities. 
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FIG. 2. (Color online) Scaling solutions for coherence length. 
The coherence length is plotted versus expansion time for (a) 
high density and (b) intermediate density. The red circles 
show results from full numerical simulations, while the solid 
blue curves show our scaling solutions using the initial con- 
densate wavefunction calculated in the trap. The coherence 
length is not well approximated with the scaling solutions 
for intermediate densities where the 0(y 2 ) expansion breaks 
down. 



mean- field phase shift at the center of the condensate, 
where it will be greatest, sets a tight upper bound. 

Therefore, an expansion valid for small deviations from 
the origin should capture most of the relevant physics. 
We expand the righthand side of (0 in a power series 
in y and keep only terms up to order y 2 . Equating the 
coefficients for the y° terms and the y 2 terms yields the 
following pair of equations for / and the A's: 
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where the a's and /3's are constants calculated from the 
initial state. 

To demonstrate the utility of these equations we com- 
pare full simulations of the three-dimensional GPE to 
approximate solutions obtained from the above equa- 
tions [23j]. We simulated the evolution of a 174 Yb (a s = 
5.6 nm) BEC released from a harmonic trap with fre- 
quencies (Q x ,Q y ,Q z ) = 2n x (50,50,20) Hz. The con- 
densate was allowed to expand for 20 ms. This is the rel- 
evant time-scale for the interferometer to measure h/m 
at the sub-ppb level described in section ITVl 

Figure [T] shows the phase at the center of the conden- 
sate as a function of time for condensates with N at = 10 4 , 
10 3 , and 10 2 (corresponding to high, intermediate, and 
low interactions). In the early stages of the expansion 
with 7V at = 10 4 , the TF result agrees well with the full 
numerical solution. This shows the initial density (peak 
density in trap of 6.2 x 10 13 cm~ 3 ) falls within the TF 
regime. However, we see that the TF result begins to 
diverge from the numerical solution after around 5ms, 
growing to a substantial deviation at the end of the ex- 
pansion. This is an important result for precision mea- 
surements. While a condensate may begin an experiment 
in the TF regime, its subsequent expansion lowers the 
density, eventually making TF no longer a good approx- 
imation. For condensates even deeper in the TF regime 
(such as the scenario discussed in section IIV|) this fea- 
ture becomes even more pronounced due to the faster 
expansion. 

For 7V at = 10 2 (peak density in trap of 6.5x 10 12 cm~ 3 ) 
the low density makes the TF approximation inaccurate 
throughout. However, the scaling solutions derived above 
maintain good validity. A common suggestion for pre- 
cision BEC interferometry is to adiabatically lower the 
in-trap density before beginning an experiment. The de- 
viation of the N a t = 10 2 results from both TF and the 
noninteracting case highlights the need to be able to ana- 
lyze interaction effects even for seemingly low initial den- 
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sities. 

The middle case of N at = 10 3 (peak density in trap 
of 2.4 x 10 13 cm~ 3 ) shows marked departure from both 
the TF and noninteracting results. The continued agree- 
ment with our scaling solution in the intermediate regime 
highlights the robustness of this technique. 

To further discuss the extent of validity of our scaling 
solutions, we now address more global properties. Den- 
sity and phase profile across the condensate generally 
show good agreement for all three interaction regimes. 
These properties individually are not usually of direct ex- 
perimental importance for atom intcrfcromctry. Instead, 
we present here the more relevant parameter of coherence 
length l c which depends sensitively on both density and 
phase profiles. For a condensate wavefunction <j) {%) we 
define the coherence length: 

fcPx<j>W<j>(x±(l c /2)x 3 ) = 1 
J d 3 x <p (%)* 4> (%) e 

For a two-arm interferometer this is directly related to 
the 1/e coherence time of the signal through the relative 
velocity of the two arms at recombination (24|. For a 
three-arm interferometer there is a similar relationship, 
up to some numerical factors related to the splitting. 

As figure [2] shows, the scaling solutions reproduce this 
global property well for high density (N at = 10 4 ) but 
show clear (25% in fig Os) deviations for intermediate 
densities (AT a t = 10 3 ). These deviations grow to 50% for 
the iV a t = 10 2 case and then decrease sharply to below 
5% for N at = 10 1 . That our scaling solutions are inad- 
equate to describe the coherence length at intermediate 
interaction strengths is not surprising given that this is 
a global quantity, and our scaling solutions are only ac- 
curate near the origin. In the small interaction case, the 
deviations fall in line with the relative size of the y 2 and 
y 4 terms in the small y expansion of the right-hand side 
of |[SJ) found using first-order perturbation theory. This 
suggests that the failure to accurately reproduce the co- 
herence length in intermediate densities directly follows 
from truncating the expansion at 0{y 2 ). We will show in 
section IIVI how our method accurately reproduces phase 
offsets for even the difficult case of a long experiment that 
moves from deep in the TF regime to weak interactions. 
Thus, extension of the scaling technique beyond 0(y 2 ) is 



beyond the scope of this work. 

III. SLOWING- VARYING ENVELOPE 
APPROXIMATION 

The class of BEC interferometers we consider typically 
contains widely disparate scales. Laser pulses used to 
manipulate condensates have durations from hundreds 
of nanoseconds to tens of microseconds [l6], while the 
entire experiment can last for tens or hundreds of mil- 
liseconds. Accurate simulation with such a separation of 
scales is computationally intensive — in some cases it may 
be prohibitively intensive. We decompose the condensate 
wavefunction into small, quasi-independent pieces, mak- 
ing simulations dramatically less computationally expen- 
sive. Our decomposition technique, known as the slowly- 
varying envelope approximation (SVEA), is well-known 
in the context of optics [25j . Its application to BECs was 
pioneered by Trippenbach et al [26j. 

The SVEA leverages one of the key experimental ad- 
vantages of a BEC: its narrow momentum spread. The 
momentum-space wavefunction of a condensate typically 
has a width well below the recoil momentum of a con- 
stituent atom due to absorption of a visible photon. 
Thus, when the condensate wavefunction is split using a 
light grating, the momentum-space wavefunction consists 
of a series of clearly separated peaks. The SVEA may be 
applied to any splitting method that creates clearly sep- 
arated momentum states, and thus is applicable to many 
possible experiments. 

We generalize the SVEA, showing that certain terms 
dropped in the usual description of the technique can be 
important at the level of accuracy needed for precision 
experiment modeling. We also identify a consistent ex- 
pansion whose small parameter gives us an estimate of 
the accuracy of this method. 

To start, we postulate a form for the condensate wave- 
function: 

4>(x,t) = e^-*-*** (10) 

where the kj's are the relevant wavevectors and the 
uij's the corresponding frequencies. We will consider 

kj = 2jk ICC where j is an integer and fc rec is the laser 
wavevector. After inserting this ansatz we reorganize the 
GPE into the following suggestive form: 



^ dt ^ 



-2^.4-3 



r 



-i(— «j— wjj +"i 2 +wi 3 )t 



„ikj -x—iujjt 



(11) 



where coordinates £j = x — Hkjt/m were chosen for the </>'s and 5 a ^b is the Kronecker delta function. The choice 
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of coordinates cancels a term of form ki 



and the 
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usual non-relativistic relation uj = hk /(2m) is used to 
cancel another pair of terms. 

Consider the Fourier transform of this equation. If 
all (pj's ("envelopes") have bounded support centered at 
zero with diameters smaller than one third [27( of the 
minimum of | fcj — kj \ , then this equation separates exactly 
into an infinite tower of equations of the form 



dt 



2m ^ 



+ 9 



^e-^^+^s-Wii-^a)* _ (12) 



h+h — li=j 

This equation differs slightly from that found in the 
literature [28]. Derivations of the SVEA refer to 
collecting "phase matched" terms into separate equa- 
tions. In the context of optics where Wj tx kj, spa- 
tial phase matching is equivalent to temporal phase 
matching. For matter waves with ojj oc k'j the 
two conditions may be considered separately. We 
have retained pieces with ki 2 + ki 3 — ki t — kj = but 
u>i 2 + lji s — U3i t — ojj ^ (i.e., we enforce spatial but not 
temporal phase matching). Such pieces can have large 
enough effects to be important in simulating precision 
experiments. Also, they are essential for including laser 
interactions in the SVEA. 

Given initial conditions satisfying the above criterion, 
this tower of equations is equivalent to the full GPE. To 
simplify we select a subset of the envelopes to keep and 
set the rest to zero. This truncation affects the envelopes 
we do keep by dropping terms from the nonlinear piece. 
We want to know how large an error this induces. 

Since we only set <f) m = if <p m (t = 0) = 0, for short 
times we can drop the kinetic term to find 
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dt h ^ 

Z2+Z3— (i=m 

For time scales short enough that the other envelopes do 
not appreciably move or expand, this can be integrated. 
The initially unpopulated terms would oscillate with fre- 
quency at least 4w roc and amplitude smaller than the 
initially populated states by the ratio of the mean-field 
energy to the recoil energy. This result gives us a good 
error estimate, since envelopes expanding and moving rel- 
ative to one another will decrease the righthand side, and 
also suggests a way to systematically improve the accu- 
racy of the approximation by retaining one or more of 
the initially unpopulated terms. 

To verify these uncertainties, we ran simulations using 
a condensate wavefunction with an initial superposition 
of three momentum states (0 and ±2Sfc rec ), all popu- 
lated with equal density. For these simulations we used 
Na condensates (a s = 2.9 nm) with N at = 10 4 and trap 
frequencies (Q x ,Q y ,Q z ) — 2n x (50,50,20) Hz. The re- 
sults from SVEA simulations were compared to results 
from full simulations of the GPE. Figure |3] shows the 
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FIG. 3. (Color online) Errors in SVEA from temporal phase 
matching. The phase at the center of the zero-momentum 
branch of the condensate in SVEA simulations is compared 
to the result of full simulations of the GPE. The blue dashed 
curve shows the SVEA with both spatial and temporal phase 
matching. The black solid curve shows the modified SVEA 
with only spatial phase matching. The initial peak of the 
dashed curve rises to 26%, but has been cropped from the 
image to make other details visible. The inset shows the solid 
curve on a scale where details are visible. 



fractional error in the phase accumulated by the zero 
momentum branch during separation. The percent er- 
ror seen in simulations with both spatial and temporal 
phase matching fluctuates with a frequency 2 x (4w rec ) 
just as would be expected for error due to neglecting a 



term 



'Wit 



-ie 



-i(4w roc )t(l 2 + l 2 ) 



The drop-off of fluctua- 
tions just after 0.2 ms corresponds to complete separation 
of the condensate branches. 

External potentials may be folded into the SVEA ap- 
proach. The potentials created by laser standing waves 
used for splitting and acceleration deserve special men- 
tion. The light-shift potential formalism shows that 
standing-wave laser pulses may be described by a po- 
tential of the form 



V(x,t) = A(x,t) 



(13) 



where A is the amplitude of the potential, which may 
depend on time and space. The use of relative detunings 
may create potentials moving relative to the lab frame. 
In the rest frame of such a potential, the e ±2zfcrcc :c pieces 
connect branches of the condensate wavefunction with 
momenta differing by ±2/ifc rcc . Since dynamics during 
a laser pulse occur at the time-scale of the recoil fre- 
quency, removing the e~ luljt phase is no longer accept- 
able. Then (fl~2"|) is modified by the reappearance of the 
term (hk rec ) 2 (f>j / (2m) , which was canceled by the time 
dependence of e~' lUjt in the original derivation. Thus, 
for modeling laser interactions it is key that temporal 
phase- matching not be enforced 29] . 
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FIG. 4. Scheme of the contrast interferometer for uu ICC and a. 
Using diffraction gratings of light, the condensate is split into 
three branches of momenta 2hk lac , 0, and —2hk lac and then 
recombined after a variable time 2T. When all three overlap, 
they create a matter-wave grating whose contrast rises and 
falls with frequency 8lu icc . 



IV. MODELING A COMPLETE EXPERIMENT 

In this section we use our theoretical tools to model 
a contrast interferometer to measure w rec and a. The 
first-generation of this experiment was reported in Q. 
This scheme makes a useful test case since it strongly 
leverages the unique advantages of BECs for precision 
measurement. 



uj Tec is extracted from the variation of (j) s i g with T. Fi- 
nally, the recoil frequency may be used to determine the 
fine structure constant a [30| . 

The first-generation experiment — which used a Na 
BEC, T — 3 ms, and momenta ±2hk ICC in branches 1 and 
3 — reported 7 ppm precision, but an inaccuracy at the 
200 ppm level [6|. Simple estimates indicated that this 
large inaccuracy likely arose from the differential mean 
field shift between the three arms of the interferometer. 
Since finite densities are needed for adequate signal-to- 
noise and since the arm imbalance cannot be arbitrarily 
controlled, residual mean- field interactions will always be 
present. This has been a strong motivation for the cur- 
rent work which is aimed at a thorough quantification of 
the mean field effect for a ppb-level measurement of the 
fine structure constant. We specifically consider use of 
Yb BECs to extend this technique to sub-ppb precision 
[3l| . Using additional Bragg pulses to give momenta of 
±40tt rec to branches 1 and 3 and T = 5 ms, the total 
phase accumulation 4> sig « 3.7 x 10 5 rad. With A at = 10 4 
the shot-noise limit is 0.01 rad. This allows accumulation 
to sub-ppb precision in less than a day. Other system- 
atic errors to this measurement can be kept within the 
sub-ppb range with reasonable choices of parameters for 
the diffracting laser beams [3l|. Within these paramet- 
ric constraints, we find that the mean-field shift may be 
reduced to below one billionth of 3.7 x 10 5 rad, such that 
atomic interactions do not preclude a sub-ppb measure- 
ment of u; rcc . 



B. Simulations of Complete Interferometers 



A. Contrast Interferometry Scheme for oj rC c and a 

In this interferometry scheme (sec figure [4} , a conden- 
sate is split into three branches with momenta 2hk rcc , 
0, and —2hk rcc by a short-duration standing-wave light 
grating (a "Kapitza-Dirac pulse") [l6|. The three 
branches evolve phase according to their kinetic energies. 
After a time T, a long-duration light grating ("Bragg 
pulse") is applied to the condensate, causing the non- 
zero momentum states to reverse direction. At time 2T 
all three states again overlap to create a matter-wave in- 
terference pattern with contrast that oscillates with time. 
This time- varying contrast is measured using a traveling- 
wave readout laser pulse with the same wavevector as the 
light gratings, which in simulations corresponds to the 
2hk ICC component of the density. 

The readout laser pulse is Bragg reflected from the 
matter-wave grating. The interferometer signal is the 
intensity of this reflection, which scales with the con- 
trast of the matter-wave grating and is proportion- 
ate to sin 2 ((</>i + <fe)/2 - cj> 2 ) = sin 2 (4u Tec t + Offset) , 
where fa is the phase of branch i (see figure |4}, 
^rcc = ^ 2 cc /(2 m ) is the recoil frequency, and ^offset is 
time-independent and does not affect the final measure- 
ment. We define (j) s ig = {4>i + </>3)/2 — 4>2- The value of 



Two full experiments were simulated to confirm the 
accuracy of our various techniques, after which we simu- 
lated the contrast interferometer of Gupta, et al [6( . The 
first allows us to test the validity of the SVEA decomposi- 
tion as compared to full numerical solutions to the GPE, 
showing that the potentially dramatic reduction of com- 
putational cost allowed by the SVEA does not diminish 
the accuracy of simulations. The second shows the effi- 
cacy of our scaling solutions in several parameter regimes 
as compared to three-dimensional SVEA simulations, for 
experiments that could not be adequately simulated on 
commodity PCs using the full GPE. 

First, we simulated a short experiment using a Na 
condensate with 1 ms of free expansion out of the trap 
((n x ,n y ,n z ) = 2n x (50,50,20) Hz) and T = 0.2 ms. 
This experiment is sufficiently short that full simulations 
of the GPE may be run in reasonable time. Figure [5] 
shows the signal from the full GPE simulation along with 
the signal from the SVEA simulation. This signal was 
generated by extracting the magnitude of the 2ftfc roc com- 
ponent of the total atomic density as a function of time. 
The phase at time 2T agrees within the granularity of 
the time steps in the simulated signal [32j. Since the fi- 
nal phase of this signal is sensitive to slight differences 
between simulations in the free propagation, condensate 
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FIG. 5. (Color online) Simulated output of a contrast in- 
terferometry experiment. The solid (black) curve shows the 
output signal simulated using the full GPE. The circles (blue) 
show the output signal simulated using the SVEA. 
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FIG. 6. (Color online) Mean-field shifts for various splitting 
parameters. The black solid line shows the scaling solution 
value for the shift in measured recoil frequency as a function 
of density imbalance for the experiment described in the text 
with £1 — 2n x 200 Hz. The blue dashed line shows the same 
but for an initial trapping frequency of Q = 2tt x 5 Hz. The 
circles show data points from SVEA simulations. The inset 
shows a closer view of the solid line. For imbalance levels of 
< 2% the mean-field shift is reduced to the lppb level. 



separation/recombination, or laser interaction periods, 
the agreement seen in figure [5] shows the power of the 
SVEA to accurately model all periods of an experiment. 
The actual experiment will consider the slope of a <p s i g 
vs. T plot in which many of these details are expected to 
cancel out. Being able to accurately model these details 
can thus improve the confidence in these cancelations, 
allow modeling of laser intensity fluctuation effects, etc. 

Several recent studies of atom-light interactions with 
varying levels of complexity [H, H3| can easily be adapted 
to describe the condensate splitting. In our framework, 
the complexity of the laser interaction model may be 
changed as needed without substantially affecting the 
models for the rest of the interferometer. For the simu- 
lations described above we used the light-shift potential 
formalism described in section Hill For the following, we 
use a much simpler model ignoring phases accumulated 
during the laser interactions, as these should not affect 
the final result of an experiment. 



The second simulated experiment is the proposed Yb 
contrast interferometer to measure a to sub-ppb preci- 
sion described above. We assume -/V at = 10 4 and the 
condensate initially confined in an isotropic trap with 
frequency fl = 2ir x 200 Hz. After a 10 ms period of 
free expansion, the BEC is split into three branches. 
The densities in branches 1, 2, and 3 are in the ratio 
(1 — x) : (1 + 2x) : (1 — x), where x is the difference in 
fraction of original density between the stationary branch 
and the accelerated branches (e.g., density splitting in 
the ratio 1:2:1 has x — 0.25). We make use of our scaling 
solutions from section |TT] to quickly probe a large region 
of parameter space to find suitable conditions for the ex- 
periment. 

To apply our scaling solutions, we break a full exper- 
iment into stages of free propagation and stages of laser 
interaction. While in section [TT] we only considered re- 
leasing the condensate from a harmonic trap at t = 0, 
the scaling solutions may be used for more general sit- 
uations by carefully choosing initial conditions. For ex- 
pansion from a trap, the initial conditions are /(0) = 0, 
Xj(0) = 1, and Aj(0) = 0. After a Kapitza-Dirac pulse, 
the condensate wavefunction has the form 



c#kd ■ 



-j3 ^exp(2i/c roc - x) + exp(— 2ik rc 



(14) 



where </>kd is the wavefunction immediately before the 
Kapitza-Dirac pulse (35j . If we ignore the interbranch 
terms from the SVEA, then the scaling equations may be 
used to continue evolving each branch of the wavefunc- 
tion forward, using O = |0kd|, /(0) = arg (4>kt>(% = 0)), 
Aj(0) = 1, and A,(0) = {2h/m)d 2 [a,rg(^ KD )]/dx 2 J \ s=0 , as 
initial conditions. 

To account for the interbranch interactions, the sim- 
plest approach is to treat each branch as though it is 
acted upon by a weak potential due to each other branch, 
ignoring the back-reaction — the change in shape of each 
branch due to its interactions with other branches. This 
approximation gives the first term of an expansion in a 
small parameter, which we describe for the case of release 
from a trap followed by some period of expansion before 
the splitting laser pulse, treating the initial expansion 
with the scaling solutions. 

If the z axis is the direction of laser prop- 
agation, we find the relative size of the ne- 
glected, second-order back-reaction effects to be 
[(Q z /oj rcc )(fj,/ (Hflz)) 3 } 1 / 2 Xy 1 for a condensate in the 
TF regime and [(O z /o; roc )(^/(ftr2 2 )) 2 ] 1 / 2 A~ 1 A- 1 in the 
small interaction limit. With this approximation we treat 
the entire post-Kapitza-Dirac pulse propagation with a 
single scaling solution and add the interbranch interac- 
tion phase at the end. These parameters are generally 
small for BEC experiments (e.g., the first-generation con- 
trast interferometry experiment @ had fl z /uj rcc < 10~ 3 ) 
showing this to be a good approximation. 

We find the fractional shift of the measured recoil fre- 
quency due to atomic interactions as a function of x (see 
figure O by calculating the phase of the signal at t = 2T 
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for runs of T = 2 ms and T — 5 ms and then finding the 
slope of 4> s ig vs. T. This agrees well with the results of 
SVEA simulations. For \x\ < 0.02 the mean-field shift 
contributes at less than the ppb level. Such precise con- 
trol of the density splitting has been demonstrated by 
Hughes et al [36l |. 

For a given atom number and available total time for 
a run, the mean-field shift is generally smaller for larger 
trap frequencies because the condensate expands much 
more rapidly after release from strong traps than weak 
traps, quickly making up for the higher initial density. 
However, as trap frequency is increased, the final mo- 
mentum spread after expansion increases, potentially di- 
minishing the advantage of a BEC's narrow momentum 
distribution. At ft = 2tt x 200 Hz, the final momentum 
spread is still less than one tenth of the recoil momentum, 
while allowing us to bring the mean-field shift comfort- 
ably below the ppb level. 

Finally, with these confirmations in hand, we simulated 
the 2002 experiment of Gupta, et al [6j, in which Na 
condensates of iV at ~ 10 5 were split with x ~ 0.25. Using 
our scaling solutions we find an expected relative shift 
due to mean-field of —3 x 10~ 4 , which compares well 
with the —2 x 10 -4 systematic error reported by Gupta, 
et al. 



V. CONCLUSIONS 

The use of a BEC source has the potential to improve 
precision atom interferometry measurements for funda- 
mental and applied research, once atomic interaction ef- 
fects are tamed. While a general rule to achieve this is to 
simply work with reduced density, doing so may compro- 
mise other aspects of the experiment such as interferom- 
eter signal-to-noise. We have presented new theoretical 
techniques to precisely quantify interaction effects allow- 
ing their analysis in high accuracy measurements at the 
sub-ppb level using BEC interferometers. 

We derived scaling solutions to the Gross-Pitaevskii 



equation valid for all interaction strengths in the cen- 
tral region of the condensate — the part which domi- 
nates most atom interferometry experiments. We demon- 
strated their validity by comparing against numerical 
simulations of the 3D GPE. 

Using our modified slowly-varying envelope approxi- 
mation, we have shown the ability to rapidly simulate 
BEC interferometry experiments that utilize standing- 
wave gratings and beam splitters. We demonstrated neg- 
ligible errors in the final read-out signal of the contrast in- 
terferometer for a when compared against numerical sim- 
ulations of the 3D GPE. Our scaling solutions adapted 
with correction terms from our SVEA analysis allow for 
precise estimation of the contrast interferometer phase. 
Large regions of parameter space (as demonstrated for 
the arm imbalance parameter) can be quickly analyzed 
to identify suitable conditions for high accuracy measure- 
ments, which would have prohibitive computational ex- 
pense using full simulations of the GPE. The mean field 
shift predicted by this method is consistent with the mag- 
nitude and sign of error reported in the first-generation 
contrast interferometer experiment [6] . 

We have shown that atomic interaction effects are not 
a road-block to a sub-ppb measurement of w roc and a 
with condensates. The presented methods can be use- 
ful in precision BEC interferometry applications where 
mean-field effects need to be well understood and the 
consequent measurement uncertainty properly accounted 
for. 
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